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Abstract 

Wc study a method to solve stationary axisymmetric vacuum Einstein equations numerically. 
As an illustration, the five-dimensional doubly spinning black rings that have two independent an- 
gular momenta are formulated in a way suitable for fully nonlinear numerical method. Expanding 
for small second angular velocity, the formulation is solved perturbatively upto second order in- 
volving the backreaction from the second spin. The obtained solutions are regular without conical 
singularity, and the physical properties are discussed with the phase diagram of the reduced en- 
tropy vs the reduced angular momenta. Possible extensions of the present approach to constructing 
the higher dimensional version of black ring and the ring with the cosmological constant are also 
discussed. 



1 Introduction 



At the level of classical theory of gravity, black holes are simple but general vacuum state of spacetime 
itself. This is a noteworthy consequence in the four-dimensional black hole physics by the uniqueness 
theorems. Nevertheless, corresponding to the fact that the number of degrees of freedom of gravity in- 
creases with spacetime dimensions, such simplicity of spacetime seems to hold only in four-dimensions, 
as was realized by the discovery of five-dimensional black rings (BRs) with horizon topology S 1 x S 2 
PQ. The BRs are in good contrast with the higher-dimensional Kerr solution (so-called Myers-Perry 
black holes), and many aspects have been studied so far including extensions of the solution, partial 
stability, role in string theory, etc [2J. 

At present, one of important aspects concerning BRs is to find out new type of "black ring" in 
higher dimensions. In more than five dimensions, we expect that not only black "rings" but also new 
types of topologies might exist. Developing a feasible effective method to construct such new black 
holes is an interesting open question. This will be also related to finding conjectured black holes with 
less symmetry than any known solutions [3]. Search for new black holes is important to understand 
the phase space of higher dimensional vacuum spacetimes. 

To tackle the issue, we need a new approach or tool. After the original investigation of BRs which 
was mainly based on guesswork (and knowledge of four-dimensional solution [1]), several systematic 
methods which can rederive the vacuum black ring solution have been developed [3 EJ El EJ EJ QUI [EH 
[T2] . An outcome is the discovery of vacuum BRs with rotation on the S 2 but without rotation along the 
S 1 circle [5]. (An independent work based on tractable coordinates is [13].) The five-dimensional BRs 
have two independent rotational planes and the two rotations can generate two independent angular 
momenta. Thus we expect doubly spinning BRs that have the rotation on both S 2 and S . However, 
the exact solution of doubly spinning BRs has not yet been discovered, and no concrete procedure to 
find it has been available. The black ring with two angular momenta is known for supersymmetric 
case, but the two angular momenta are not independent conserved chargespH], 

The present paper aims at exploring a new approach to search for new type of black holes. A 
numerical approach has been employed and developed to discover static axisymmetric solutions over 
the past several years [151 EE El 021 02] • The advantage of numerical approach is that it will be widely 
applicable once the basic scheme and techniques are established, and therefore it will be interesting 
to extend the previous approach applied for static axisymmetric spacetimes to stationary ones. To 
construct a solution numerically, it needs adequate understandings of expected solution. However, as 
we will see, BRs are rather complicated in numerics to hold a consistency over whole system, and 
boundary conditions for unknown solutions are not necessarily clear a priori. In numerics, we need 
to fix all necessary information before starting the calculation. This is in contrast with an analytic 
approach, in which we can impose necessary conditions after we find a general solution, roughly 
speaking. 

In this paper, we adopt the canonical form of metric for stationary axisymmetric spacetimes 
[201 [T2] as suitable coordinates for numerics, and we will work on five dimensions. Although the so- 
called ring-coordinates [2] is the most adopted coordinate system for BRs, it is not rectangular at the 
asymptotic infinity (boundary) and less useful. As a good example to test such approach, we focus on 
the issue of doubly spinning BRs. We formulate the doubly spinning BRs in a numerically manageable 
manner, and solve it in a perturbative way. Based on the solutions, some of physical properties of 
the doubly spinning BRs are discussed. This specific example helps us to clarify subtle issues related 
to numerics, and we discuss how we can deal with fully nonlinear stationary solutions in the present 
approach. Based on these, we also discuss a possible extension to find higher dimensional BRs with 
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new topologies and BRs in AdS. 

The paper is organized as follows. In the next section, we introduce the canonical form of metric 
and explain necessary ingredients. In Sec. [3l we formulate doubly spinning black rings in a generic way 
and provide the rod-structure which works as a boundary condition. In Sec. [U we solve the doubly 
spinning BRs in a perturbative way and discuss the physical properties in details. Several possible 
extension of the present approach are discussed in Sec. [5j The final section is devoted to a summary. 



2 Stationary axisymmetric spacetimes 
2.1 Canonical form of metric 

Following to the formulation in [2D] (and [21] for static case), we introduce the canonical form of 
metric and collect necessary ingredients for later discussions. Stationary axisymmetric D-dimensional 
spacetimes can be put in the following form, 

D-2 

ds 2 = 9abdx a dx b + e 2v {dr 2 + dz 2 ). (1) 

a, 6=1 

If this spacetimes has D — 2 commuting Killing vectors, the metric is in a special class. For such cases, 
the canonical coordinates in which the coordinate r satisfies 

r 2 = det(-<?), (2) 

are available. This canonical form of metric simplifies the Einstein equations considerably. The 
vacuum Einstein's equations are 

2e 2u R ab = -V 2 g ab + g dc (d r g ad d r g bc )+g dc (d z g ad d z g bc ), (3) 
Rzz = -V 2 v-\g ab g cd d z g a cd z g bd , (4) 

Rrz = -d z v - \g ah g cd d z g ac d r g bd , (5) 
r 4 

12 1 

R rr - Rzz = -n + -d r v + -g ab g cd (d z g ac d z g bd - d r g ac d r g bd ), (6) 

where we have introduced the operator V 2 = d 2 + d 2 + ^d r , which works as the Laplace operator 
in three-dimensional cylindrical coordinates. The first two differential matrix equations with the 
constraint ([2]) are the basic field equations to be solved. The last two equations are the constraint 
equations. It can be shown that the constraint equations are satisfied as long as ([3|) and (jl]) are 
satisfied. If we find a solution of ([3]), v(r, z) will be immediately obtained from and thereby a 
complete solution of the Einstein equations will be obtained. 

Putting x 1 = t, x 2 = <fi, x 3 = ip where both x 2 and x 3 are periodic coordinates, the five-dimensional 
Minkowski spacetimes is 

2 

ds 2 = -dt 2 + (Vr 2 + z 2 - z)d(p 2 + (yV 2 + z 2 + z)d^ 2 + - 6 (dr 2 + dz 2 ). (7) 

2Vr 2 + z 2 

The constant parameter e is related to the periodicity of 4> and ip, and they are given by A(j> = Aip = 
2ire. Appropriate rescalings of r,z and the angular coordinates allows us to take e = 1, and it is the 
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Figure 1: Canonical coordinates for Minkowski spacetimes ([7]) and black rings ([8]). 

standard form of the Minkowski spacetimes without conical singularity. However, the above form of 
Minkowski metric will be useful when we consider the spacetimes of black rings in which the periodicity 
is not necessary to be A(p = Aip = 2ir. The coordinate system is illustrated in Fig. [TJ 

In five-dimensions, the exact solution of rotating BRs can be expressed in the canonical form of 
metric. The spinning black ring is given by 

W. r 2 9 2 H (R 3 + z- k 2 )(R 2 - z + ck 2 ) 

gtt — 9H — 1 ) 9ipip 
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W + gtt9ij>ij> 9tt Ri-z- ck 

AcKy/QTcj R 3 -R! + (l + c)k 2 2u _ (1 - c)Ri + (1 + c)R 2 + 2cR 3 

9t<t> v 7 !^ W+ 6 + 8(1 - c^RrfzRz ' { ) 

and = = 0, where we have introduced 



Ri = \Jr 2 + (z + ck 2 ) 2 , R 2 = \Jr 2 + [z - ck 2 ) 2 , R 3 = yjr 2 + (z - n 2 ) 2 , 
W± = (1 + c)Ri + (1 - c)R 2 - 2cR 3 ± 4ck 2 , 

Here, < c < 1 and k > are parameters. Roughly speaking, the parameter k sets the scale for the 
solution so that it measures the radius of the ring, c measures the ratio between the horizon radius of 
S 2 and the radius of the ring. So smaller values of c corresponds to thinner ring. In the limit of c — > 
and k — ► 0, the Minkowski metric is recovered in the form ([7]). The BRs have no conical singularity if 
the periodicity of angles takes 

27T 

A<f> = Ail> = 2vre (f)9) = (9) 

Vl + c 2 



2.2 Rod-structure 

The stationary axisymmetric solutions expressed in the canonical form of metric can be uniquely 
characterized by its boundary condition on the z-axis, known as rod-structure. The rod-structure 
can be rephrased as physical requirements/conditions of a black hole solution, so that it plays an 
important role to find a new solution. 
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Let us briefly discuss what the rod-structure is. The z-axis can be divided into the n + 1 intervals, 
[—00, zi], [zi, Z2], ■ • • , [Zn, doc]. A necessary condition for a regular solution is that precisely one eigen- 
value of 5 a fe(0, z) becomes zero for a given z, except at the isolated points between the intervals. The 
n + 1 intervals are called rods of the solution g ab , and the rod-structure for a solution is defined by 

\\mg ab (r,z k )v h , k s = 0, (10) 

with non-vanishing vectors vha (k = 1, • • • , n), on the fc-th interval (rod). Each vector associated 
to a rod has its own direction. If \\m r ^Q g a bV a v b /r 2 is negative (positive), the direction of the rod is 
timelike (spacelike). A timelike rod corresponds to a horizon, while a spacelike rod corresponds to 
a compact direction. Moreover, a semi-infinite spacelike rod corresponds to an axis of rotation, and 
the associated coordinate is a rotation angle. These properties of rod and direction provide useful 
information of spacetimes. 

Taking a closer look at the EOMs near r = 0, we find that a solution g ab behaves as g a b = 
0(1) + 0(r 2 ) + • • • , where the coefficient of each term will be a function of z. Except the isolated 
points between the rods, g a b stay finite at r = in general, but some metric components vanishes 
according to the rod structure. Suppose gu behaves as 511 « n(z)r 2 on some rod. Then, from R rz = 0, 
we find 

e 2 " = j 2 n(z) + 0(r 2 ) (11) 

to leading order for r = 0, where 7 is a constant. Then the relevant part of metric is 

ds 2 ~ n{z)[r 2 {dx 1 ) 2 + 7 2 dr 2 ]. (12) 

If x 1 is a spacelike coordinate, a conical singularity appears on the axis, unless x 1 is a periodic 
coordinate with period 

Ax 1 = 27T 7 . (13) 

Therefore, for a rod in a spacelike direction, a specific periodicity is required to avoid a conical 
singularity. The same argument also holds for a timelike coordinate. If x 1 is a timelike coordinate, 
there is a horizon at r = 0. By a Wick rotation, the Euclidean time must be periodic and it means that 
the temperature T associated to the horizon is T = l/2ir^. For stationary case, however, the direction 
v b will have multiple non-vanishing components and we have to take into account it. Furthermore, 
the BRs have no Euclidean counterpart that is real without conical singularity. 

3 Formulation of doubly spinning black rings 

The canonical coordinates have rectangular boundary, and they are suited for numerics. Under the 
canonical form of metric, the field equations are non-linear coupled equations with elliptic second-order 
partial differential operator. Besides, the boundary conditions are well imposed by the rod-structure, 
once we fix the type of black holes. Thus we will be able to apply a numerical method for stationary 
cases. In this section, we formulate the black ring with two spins and discuss the necessary boundary 
conditions, including the basic strategy to solve the equations. 
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To deal with general black rings carrying two spins, we define the general metric in the following 
form; 

e = e M u , g H = g^'M^, = g^ J M^, 

m = gf^ Mt^M^, g tl p = gf^Mt^M^, g H = M^. (14) 

Although we define g tt = g^M tt , it is solved from the algebraic relation (JSJ). gf^ is defined by 

(bg) _ ( 2CK 2 V {hg) 

Here, M a b(r, z) and M u (r,z) are unknown functions which we solve. We adopt the Emperan-Reall's 
solution ([8]) as the "background" metric g^J' , so that the Emparan-Reall's solution is recovered in the 
limit of 

Mttj, , , Mfap , M v -» 1, 

M^M^-tO. (16) 

Our metric ansatz, in particular and gt^, are chosen in order to yield simple expression for physical 
quantities, such as the angular velocity of horizon. Besides, as we will see below, an advantage of 
adopting the above ansatz is that the full metric will automatically satisfy the desired rod-structured 
as long as M a {, are smooth functions. 

From Eqs. ([3]) and ([4]), the coupled non- linear PDEs for the functions M a b are symbolically written 

as 

V 2 M = S(M), (17) 

with the non-linear source term S, where M = M a ^ or M v . To solve this system as a boundary 
value problem, we specify the boundary conditions by explicitly analyzing the rod structure and the 
asymptotic behavior for the assumed general metric §14 



3.1 Asymptotic analysis 

Let us first specify the asymptotic behavior of the solution. In the following new coordinates, 

r = p 2 sin#, z + ck 2 = p 2 cos 9, (18) 

the metric asymptotes to the five-dimensional Minkowski space ([7]) in the limit of p — ► oo. Assuming 
the general periodicity of A</> = Aip = 2ire, the subleasing order terms are found to be 

4M . 4(M + 7 ? )sin 2 (#/2) . 4(M - rj) cos 2 (6/2) 

^""W 3^ ' 3^ ' 

^sin^/2) 2J,,co S 2 (e/2) 5a M ~- C -^ (19) 

Here, M and J are mass and angular momentum, respectively, and £ and ij are constants. The 
periodicity e is related to the leading order of e 2L> , which is given by 
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6 -v +0(1) ' (20) 
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From these behaviors and those of background metric ([8]) , the metric variables M ao must obey the 
following asymptotic decays, 



M ab ~ 1 + 



Oiab 



Mab ~ a ah + 0(p~ 
sin 2 9 

M^ ~ a^j, 2 , 



I, (M t4> ,M til ,M„ 



(21) 



were a a h are constant coefficients. These asymptotic behaviors are necessary to be imposed as the 
Neumann boundary conditions at the outer boundary. From each coefficient, the physical quantities 
can be read off as 



Af = — — a tt , 



M 



3vr^ 



4ck 2 
(1-c) 



+ a^V 



2vre 3 K 3 c(l + c) 3 / 2 3 2 4 

= — , J2 = Tre c k a^, 

where att yields att 4ck 2 /(1 — c) in the limit of J2 — > 0. 



2 2 

" £ (6ff)"" 



(22) 



3.2 Boundary conditions and Rod structure 

We analyze the boundary conditions of the general BR metric and specify the rod-structure. Since the 
geometry of the general BRs with two spins is the same as the black ring with single spin, the basic 
rod-structure should be also the same. The rod-structure is classified into four parts: the semi-infinite 
spacelike rod [— 00,— ck 2 ], finite timelike rod [— ck 2 ,ck 2 ], finite spacelike rod [+ck 2 ,k 2 ], and the semi- 
infinite spacelike rod [k 2 , 00]. Near the r = 0, the metric components behave as g a b = 0(1) + 0(r 2 ), 
so that the basic boundary conditions are the Neumann boundary condition: 



d r g a b = 0. (r = 0) 



(23) 



In addition to this condition, some types of Dirichlet boundary conditions are necessary. An important 
point is that the boundary conditions, which can be translated into the rod-structure, is determined 
by requiring physical (or regularity) conditions on the z-axes. 



The semi-infinite spacelike rod [k 2 , 00] corresponds to the 0-axis. This means that 



9t<f> — 9<f>ip 







(24) 



should be satisfied there. This physical requirement is equivalent to having a rod in the direction 
v = (0,1,0), i.e., the d/dx 2 direction. From the property of background metric, and g tt p 
always satisfy the above condition. The non-trivial condition is 



(jnj) 



0. 



(25) 



As discussed in Sec. 12.21 the periodicity of angular coordinate is given by A(f> = 2tt lim r ^o \/r 2 e 2u 
which yields 



2tt 



My 



VT+^y M,, 



-axis) 



(26) 
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Since the periodicity does not change along the axis, the above result means M v jM^ = const, 
on the axis. In the asymptotic limit (p — > 0), we obtain — > 1 and M u — > const., and so the 
asymptotic value of M v determines the periodicity (see Eq. ([20]) ). The condition (f26j) can be 
used as a Dirichlet condition for M u , once its asymptotic value is fixed. 

Similarly, the semi-infinite spacelike rod [ck 2 ,k 2 ] corresponds to the axis of ip which locates 
inside the ring. We call this portion of axis the inner axis. In this case, the physical requirement 
(and the rod structure) gives 

9iptp = dtip = 9<t>i> = 0, (27) 

and we find = 0. The direction of the rod is v = (0, 0, 1). Note that the background metric 

dty* ls cnosen so as to satisfy this condition automatically. With the same argument above, the 
periodicity of this axis is found to be 



27T / M u 

Ai> = ^===J—^, (V-axis) (28) 

v l + c z y 

to avoid a conical singularity. This condition can be also used as a Dirichlet boundary condition 
for M v . 

The finite timelike rod can be put at [— ck 2 ,ck 2 ], and it is the horizon. From (|10p with v being 
the vector v = (1, u;^, aty), we find that the angular velocities on the horizon are 



■yJ\ — C 

w = 7. n M t<P> = M H- ( 29 ) 

lK\J 1 + C 

Here we have required 

= 0. (30) 

The physical angular velocities are given by 

2-7T 2-7T 

As we notice, Q,^ and fi^ are constant over the horizon as long as and M^j, are constant. We 
will be able to use these constants as parameters of specifying a solution. This is the reason why 
we have adopted the metric ansatz in the form of (|14p . Note that physically the rod-structure on 
the horizon is equivalent to finding a null Killing vector £ a = (dt) a + u^(<9</,) a + uj^(d^) a , where 
lo^ and uty are defined by io b = -g t b/9bb with = 0. 

From the generator of the Killing horizon, the temperature is 

(32) 



2-7T7 ' 



7 = 



(l-c)T(TT?) 
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Here, I/7 is the surface gravity. Prom the constancy of surface gravity on the horizon, M^M^M v \ r= § 
is constant over the horizon. Then the entropy is given by 

2c 2 k 3 A0Ai/> 



S = Z-Z- s/M^M^M. (33) 

(1 - c)vl + c l 

Note that the entropy and temperature follow the Smarr relation and the first law with other 
physical quantities: 

-M = TS + V^iJj, dM = TdS + J2 n i dJ *- ( 34 ) 
3 

The infinite spacelike rod [—00, —ck?] becomes the axis of tp which locates outside the ring. The 
direction of the rod is v = (0,0, 1), which means 

9ip4> = 9tip = 9<pip = 0. (35) 

Again, we obtain a non-trivial condition = 0. To avoid a conical singularity on this axis, ijj 
needs to have the period 



The above analysis of rod-structure gives us additional information on the boundary conditions, 
and all of these conditions are sufficient to solve the Einstein equation as a boundary value problem. 
Our choice of the background metric is chose so that the full metric automatically satisfy the 
Dirichlet boundary conditions (|24p . (|27p . (|35p on the axis except for g^ , in addition to providing simple 
expression for the rotational velocity on the horizon. The difference in the rod-structure between the 
BR with single spin and the doubly spinning BR is the direction of the vector v on the horizon. 



4 Perturbative approach 
4.1 First order perturbation 

We are ready to solve the field equations as a well-defined problem. Although it is very challenging 
to develop an elaborate numerical code that solves the nonlinear equation (|17p directly, we discuss a 
method to solve the system perturbatively. Passing this approach, we can confirm the validity of the 
formulation and the specification of the rod-structure discussed so far. 

Based on the behavior (|16|) near the background metric, the metric functions are expanded as 

00 

M ab = m^+J2^, (37) 

n=l 

where uj^ is an expansion parameter and the first term is Tn& = or 1. For instance, 

M H = 1 + ti^rn^l H , 

M H = w^mJJ + • • • . 



(38) 
(39) 
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For brevity we will omit hereafter the superscript [n] denoting the order of perturbations whenever it 
is obvious. From (|29p (on the horizon), the expansion parameter tly, is just the second angular velocity 
on the horizon at the leading order, 

^ip ~ £y>) (40) 

and we do not need to distinguish them hereafter. This is the most natural expansion parameter in 
the present issue because the angular velocity vanishes in the limit of black rings with single spin and 
will remain small finite value even for maximum rotation. 

At the first order of uj^, the coupled equations are symbolically written as 



V 2 m^ = 0, (42) 



3> t1p = 0, (41) 



X7 2 m u = Sim^mt^jTn^j,). (44) 

Here S is a source term depending on variables m a b- 7 and Q represent a general second-order elliptic 
operator, e.g., 

9 = Adl + Bd r + Cd 2 z + Dd z + E, (45) 

where the coefficients A,B, - - are 2x2 matrices of functions of {r, z}, which contain background 
metric. It is remarkable that the equations for {m^, m^w,} decuples from the others. Thus by solving 
(|41|) . we can introduce the second angular momentum into the BRs with single spin. Eq. (|42f) is 
solved analytically, and one finds that all non-trivial solutions do not satisfy the asymptotic boundary 
condition. Hence, we have m^p = 0, and Eq. (H3j) is also decoupled from the others. 

The equation (|43p with = has a regular solution. The angular velocity oj^ oc on the 
horizon is a constant parameter, and once it is provided, the solution is uniquely determined. In fact, 
if we set = on the horizon, only the trivial solution is found numerically. Note that the operator 
Q are invariant under a constant shift of — ► + const., m%& — > mts — const. This constant 
mode will be fixed by the asymptotic condition — ► 0. 

The meaning of the independent perturbations of (|43h is clearly that the mass and angular momen- 
tum of the background metric can be changed infinitesimally irrespective of the second spin. This is 
consistent with the continuous sequence of BRs with single spin. We are not interested in these modes 
so that we turn off these perturbations. In this case, m v has no regular solution except a constant 
mode, and so we also kill this perturbation. After all, we take 

rn [ l\=m [ l\ = m [ ll=m^=Q, (46) 

at the linear order. 

Even with these simplifications, the above PDEs for the perturbations do have complicated coef- 
ficients which are functions of {r, z}, and it is impossible to carry out separation of variables, except 
for m« and m v . We solve the PDEs by employing a numerical method. An example of numerical 
solutions is shown in Fig. [2] for a particular set of c and k. From a sequence of numerical solutions, 
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Figure 2: Numerical solution of and in the coordinates (fl~8|) for c = 0.18 and k = 2.9. The 
thick and dashed lines on the p-axis (6 = 0) represent the horizon and </>-axis in the {r, z} coordinates, 
and the portion between them corresponds to the inner axis. The outer axis is at 6 = w. The 
asymptotic boundary is at p « 12, which is not covered in the above figures. 



we read off the angular momentum J^. In Fig. [3j we show how the second perturbative angular mo- 
mentum changes with respect to the background parameter c. To eliminate the expansion parameter 
aty and k, we combine the following two quantities at the leading order, 

J 4c(l + c) 2 ttk 4 



(C 3 - C 2 + C - 1) 2 ' 
Jw, 7TC 2 K 4 

?V = (TT^F^- (47) 



Figure [3] shows that the second angular momentum increases with respect the background angular 
momentum J a, as c increased until c < 1/2 (thin black rings). Notice that the background angular 
momentum J<p (and J^/fl^) increases monotonically as c — > 1 (limit of thick black rings). Then the 
above result means that for a fixed angular velocity fty the perturbative contribution of ■/■/, becomes 
relatively larger as c increased. For the branch of fat BRs (c > 1/2), J^/Vt^ diverges in the limit of 
c — > 1, and so the relative contribution of J^/fty decreases as long as J^/Q,^ is finite. 

4.2 Second order perturbation 

It is interesting to consider the backreaction from the second angular momentum on the background 
BRs. It will allow us to check the consistency of perturbations and will give useful information about 
the doubly spinning BRs. The relevant Einstein equations at the second order are 

v2m !w> = s ( m H>, mw), (48) 

/ I 2 ! \ 

m 



Q [2] I = S(m^j,,mt^,m^i,), (49) 
V m t<t> ) 

V 2 m[? = S(m^, m t( i„m^,mtii>,m^), (50) 
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Figure 3: Plot of (J^A^X-V^)" 



where S represents source terms which contain quadratic of linear perturbations and /or second order 
variables. The inhomogeneous solutions of these equations carry the backreaction of second angular 
momentum. These equations will be solved numerically in descending order. 

There are subtle issues in (numerically) solving these equations. Let us first discuss the equation 
(|49|) . The Neumann condition d r m a i, = are the basic boundary conditions for and mw, at r = 0. 
The exception is on the horizon, which is constant over the horizon. Since and have 
the homogeneous solutions, the constant value on the horizon can be set arbitrary by adding the 
homogeneous solution. However, this ambiguity is removed by noticing that the periodicity on the 
0-axis and that on the ip-&xis should be the same, to achieve the balance between forces in the ring. 
Expanding Eqs. (J26]) and p8j) . 



Aip 



2tt 



2vr 



V / TT? 



1 + -^-(m v - m^) 



Assuming the continuity of m„ at the boundary between </>- and ■z/'-axes, we obtain 

= mw (r = 0,z = k 2 ) 



(51) 



(52) 



This condition fix the overall amplitude of m^; the homogeneous solution for (|49p is superposed on a 
specific solution of (|49p so as to satisfy the above condition. With this condition for m^, the angular 
velocity ua, oc ra^ is determined automatically. 

Next, we discuss about m u . Let us recall that this variable is defined by 



Jlv 2v, t 



e*«M v = e zv ^) (l + uj^m v + ■■■), 
and its asymptotic value is related to the periodicity parameter e, 

m v = fi„ + 0(p~ 2 ), 

2 2 n , 2 ~ \ 
6 = e {bg)\ l +U i ,U v ), 



(53) 



(54) 



where a u is defined by a u = 1 + w 2 a?;,. Although the solution for m u is easily found by inverting the 
operator V 2 , we encounter a technical deficit in the present approach whenever we adopt numerical 
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solutions. The issue is how we fix a constant mode involved in m„, which directly affects a u . We 
notice that it is free to add a constant to a particular solution of m v , which corresponds to changing 
the amplitude of M u , or adding a constant in is, 

v — > v + const. (55) 

Adding a constant mode is completely consistent with all the boundary conditions for e 2u , and no 
boundary conditions (and the PDE) can fix it. The value of m u (or u) is related to the periodicity of 
axes, which enters into the physical quantities, and so it has a physical meaning. However, as long as 
we solve the PDE numerically, we can add an arbitrary constant to each obtained solution. One might 
think that we are able to fix the ambiguity through the Smarr relation between physical quantities. 
But this does not work because the periodicity e is factored out as an overall factor, and the Smarr 
relation does not care about the periodicity [22]. 

What is the origin of this problem? Suppose that we find a full analytic solution for the PDEs, and 
consider how we fix the amplitude of e 2v ' . The amplitude will be fixed by requiring that the periodicity 
of angles becomes, e.g., 2tt in the (continuous) limit of Minkowski spacetimes. This is the only way we 
can actually do for the black rings with single spin. On the other hand, however, if we find solutions 
numerically, they constitute a discrete sequence of solutions and each numerical solution can carry its 
own (arbitrary) periodicity of A<p = Atp. Since the periodicity (or arbitrary added constants) has no 
relation with each other, we have no principle to fix it a priori. 

Except for this problem, there is no other ambiguity in solving m v . To fix the constant mode, we 
will set a v = when we solve m v . Once we fix the constant mode, the boundary condition on r = 
are just Dirichlet conditions derived from (126f) . (|28f) . (I32j) and f)36f) . where we use the fact that the 
temperature and periodicity are constant. For instance, Eqs. (f26|) and (f28|) gives 

m" =0 = m^, (0-axis) 

m"~° = m^j,., (if)— axis) (56) 

where m° =0 is a solution with a v = 0. 

Based on these considerations, we have solved the second order perturbations numerically. An 
example of solutions is shown in Fig. [U The perturbative solutions are constructed with respect to 
the thin black rings (c < 1/2) and their physical properties are discussed in the next subsection. For 
the fat black rings (c > 1/2), the control of the homogeneous modes in Eq. (I49h becomes harder, and 
we could not construct a solution on this branch. Thus we restrict ourselves to the doubly spinning 
thin BRs in the next section. From Fig. we notice that the solution m v has a sharp peak at 
the boundary between the horizon and inner axis (at 6 = 0). This non-smooth behavior would be 
inevitable, because the background metric itself is not sufficiently smooth in the canonical coordinates; 
e 2 "( b 9) diverges at z = ±ck 2 , k 2 on r = 0. 

Before ending this section, we give comments on the accuracy of numerical solutions. The Smarr 
relation (]34j) is expanding upto second order and no linear order term appears. As discussed above, 
the Smarr relation does not suffer from the issue of the constant mode, this relation is useful to check 
the consistency. We have confirmed that the Smarr relation at the second order is well satisfied within 
2% relative errors. Furthermore, we also found that the relation still holds within the same error even 
if we do not impose the condition (|52p . This is consistent with the fact that the Smarr relation is 
satisfied even if the force balance condition A(j) = Aip is not enforced as discussed above. 
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Figure 4: Numerical solutions of m^, m^, mw, and m^, corresponding to the solutions at linear 
order in Fig. [2j 



4.3 Perturbative doubly spinning black rings 

By the perturbations upto second order, we have incorporated the backreaction from the second spin. 
The backreaction affects the spacetime geometry, and so we can investigate the physical properties of 
the doubly spinning black rings based on our solutions for thin black rings. 

All physical quantities depend on the periodicity e, so that they are not free from the ambiguity 
of a v (or the constant mode in m u ). To avoid the ambiguity, we adopt the so-called reduced spin 
variable j and reduced horizon entropy, 



■2 



J 



27tt J% + Jl 

32 M 3 ' 
3^3 S 

Here we have introduced the reduced spin by a sum of squared angular momentum. This is useful 
variable to draw a phase diagram of doubly spinning BRs in a two-dimensional plane, instead of in 
a three-dimensional space, j 2 does not depend on the periodicity e by its construction, while the 
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reduced entropy does depend on it. However, expanding these quantities upto second order, we find 



■2 -0 2 / -2 -2 \ 

3 = gc +< ^Cfy[2]+W' 

a = 2v/c(l-c)+u;j5 [2] , (58) 
where the second order corrections are 

. 2 «^^ 2 (1 - C) 3 

%2] - 32 
(1 + c) 3 

4[2] = -3-2-^2" (8ck 2 «^ - 3(1 - c)(q^ + a^)) , 

f7F\ 7\(m Mm + ™«=o 3(l-c)(a^ + q^) ^ 

s [2] = V c \ l ~ c ) I m VV> + m 00 + m u -^2. J> 

and we also use [j 2 )[2] = i|r 2 i + J^pi- j|[2] ana ^ ^[2] represent contributions from respective angular 
momenta. As it shows, the dependence of a v in s is higher order in perturbations so that it is free 
from the ambiguity at the leading order 0(cj3,). For the MP black holes, the reduced entropy is 



smp = 2yi - 4 - 4 + ([1 - (u - u) 2 ] [1 - (u + u) 2 ]) 1/2 > (59) 

and the horizon exists only when j 2 + j| + 2\j^j^\ < 1. 

In Fig. El sp] and (j 2 )[ 2 ] are shown. Both of them are negative, implying that the doubly spinning 
black rings have less entropy for given mass. It is a general process to decrease the horizon area as a 
spin increases, as are the cases in the Kerr and MP black holes. However, we have no reason to expect 
smaller reduced angular momentum, and so a question is why (j 2 )[2] is negative. In the evaluation 
°f (j 2 )[2]i &2\ gives small positive contribution while gives dominant negative contributions, i.e. 
|j|j 2 j| 3> &2]- Thus the negative value of (j 2 )[2] comes from the backreaction to j| and we have 
(j 2 )[2] ~ 3a, < 0- By comparing with the MP black holes (Fig. [5]), the present result might suggest 
that there is a maximum reduced angular momentum j 2 and its maximum value is achieved by the 
single spinning BRs. This does not necessarily mean that the maximum total angular momentum 
J 2 = J 2 + J 2 is saturated by the single spinning black rings. We expect that the angular momentum 
of the single spinning black rings provides lower limit of the total angular momentum, because the 
second angular momentum would have no significant role for the radial balance of forces in such black 
rings. Therefore, the negative reduced angular momentum means that the increase of mass due to the 
backreaction is more rapid than the increase of spins in these perturbations. 

Finally, we discuss the ergosurface of the doubly spinning BRs. The ergosurface at which d/dt 
changes from timelike to spacelike is shown in fig. [6] for the background BR and the doubly spinning 
BRs. The doubly spinning black rings have larger ergosurface, and the surface tends to spread outside. 
This behavior is opposite from that observed for a black ring with a rotating 2-sphere [7] which has a 
conical singularity. 

Related to this geometrical property, we give a comment on closed timelike curves. The black ring 
with a rotating 2-sphere [7] has closed timelike curves, i.e. g^p < 0, in some parameter space. On the 
other hand, it is clear that closed timelike curves are not allowed as long as we work on the perturbative 
expansion ([37|) . since g^ and g^ cannot change the sign by perturbations. Furthermore, the positive 
amplitudes of and means that the actual perturbations has no symptom of closed timelike 
curves. 
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Figure 5: (a) The second order correction in s, s/j 2 and 10" 

for the black rings (solid) and Myers-Perry black holes (dot-dashed). The curves for MP black holes 
correspond to J^, = 0, J$ = J^/4:, and J$ = cases. Since the existence of horizon requires 
j| + + 'AiipjA — 1> the latter branch ends at some point. Each dashed lines emanating from the 
sequence of thin black ring indicates the sequence of doubly spinning black rings. 



5 Extensions of the approach 

So far, we have focused on the perturbations to understand how we solve the boundary value problem 
(|17p for a given rod-structure. Besides, by the example of doubly spinning black rings, we have made 
it clear that the rod-structure can be actually specified by requiring physical properties of an expected 
black hole. 

Clearly, it is possible to extend the perturbative solutions to fully nonlinear solutions by directly 
solving Eq. (|17p . Although the fundamental problem of the constant mode in m u remains, all boundary 
conditions explained in Sec. 13. II and 13.21 can be applied, and an iterative scheme to solve elliptic PDEs 
are available. These techniques are also useful to discuss other interesting systems. For instance, 
several static axisymmetric solutions (black holes or bubbles) will be extended to stationary ones 
[211 123] . or as an extremely exotic possibility, it might be possible to construct a black hole (or black 
ring) put on the center of the Emparan-Reall's black ring. 

Another interesting possibility to apply fully nonlinear numerical method is the search for black 
rings in higher dimensions. There have been no known feasible approaches to find such objects so 
far. Here we discuss a possible approach to construct such object based on numerics. A crucial point 
to construct higher dimensional black rings and other objects is to find an appropriate coordinate 
system, as is often the case in general relativity. In particular, it is best in numerics to find rectangular 
coordinates in which all boundaries coincide with the coordinate lines. For higher dimensional black 
rings, the canonical coordinates ([2]) is no longer available in D > 6, but the general conformal form of 
the metric is still valid. 
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Figure 6: Ergosurfaces for the background BR with single spin (inner solid curve) and the doubly 
spinning BR (outer dashed curve with oj^ = 0.3). To compare the two ergosurfaces, we draw the thin 
dotted line by just rescaling the inner solid curve for the background black ring. The doubly spinning 
black rings have larger ergosurface, and the surface tends to spread outside. The thick line on r = 
represents the horizon. 



The metric ([T]) can be written as 

D-n-2 

ds 2 = 9ap(r,z)dx a dx p + e 2X ^daf n) + e 2u ^ z) {dr 2 + dz 2 ). (60) 

a,/3=l 

Here da 2 n ^ = r jijdx l dx'-> is the metric of n-sphere. This form of metric will be suited for describing a 
black object with S l x S- 7 topology. In fact, the Minkowski spacetimes for D > 5 can be embedded in 
this form: 



9apdx a dxP = -dt 2 + {Vr 2 + z 2 - z)dQ 2 D „ 

a,/3=l 



e 2X = yV 2 + z 2 + z, (61) 

-,2u 1 



e 



2Vr 2 + z 2 ' 



The corresponding ring coordinates were first discussed in [23]. For D = 4, the metric is the so-called 
Papapetrou form, 

ds 2 = -dt 2 + r 2 d(f> + dr 2 + dz 2 . (62) 
For D = 5, making the coordinate transformation 

r = ift 2 sin 26, z = ^R 2 cos 26, (63) 

with 0<6<7r/2, we get the Minkowski metric in spheroidal coordinates. 

ds 2 = -dt 2 + Jl 2 cos e 2 dip + d'R 2 + Jl 2 (de 2 + sin 2 6# 2 ). (64) 
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For D > 6, the same coordinate transformation as (|63p is available. 

Taking these Minkowski metric as a background metric and assuming a general form of metric, 
like Eq. (|14p . the PDEs are reduced to coupled elliptic PDEs with the Laplace operator V 2 in 
three-dimensional cylindrical coordinates. The canonical condition ([2]) is not available anymore, and 
important properties of rod-structure is also unavailable. For instance, as we notice from the above 
example, the principle of "one eigenvalue of g a b becomes zero for a given z" is broken explicitly on 
axes, since multiple angular components, say <if2 2 , vanishes at the same time. However, we can still 
require physical (or regularity) conditions on the axes and other portion of r = 0. Then it is possible to 
set up a problem of constructing higher dimensional black rings as a well-defined (numerical) problem. 
(This is analogous to finding a black hole localized on the KK circle in D = 5, whose construction in 
numerics had been done without relying on the rod-structure |2 HI16^ [T7]). Subtle issues will come with 
rotation of the rings. For instance, it will be necessary to factor out angular variables related to sphere 
in the off-diagonal metric components, and the associated spheres may not be round sphere. Besides, 
we may also worry about the ergosurface, at which vanishing of gu may cause numerical problem. 
Nevertheless, higher dimensional BRs are great challenging issues in higher dimensional gravity, and 
further research is necessary to tackle the issues. 

Finally, we end this section with another interesting possible extension of the present approach. 
From the point of view of numerics, the cosmological constant is not an essential problem once a good 
coordinate system is provided. Here we provide the metric of five-dimensional AdS in the conformal 
form, 



X 



6 ' y x 2 

A 



r 2 + z 2 - zj dcj) 2 + (\A~ 2 + z 2 + zj dijj 2 + e 2u (dr 2 + dz 2 ) 



1 + — x 2 



X~ = 2yr 2 + z 2 . (65) 
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where e 2l> is the same as in (|6ip and A < is the cosmological constant. The domain of x 1S 
X = [0, \/24/|A|], and so h = [0, oo]. This metric is found by coordinate transformation from the 
static coordinates of AdS. 



6 Summary 

We have discussed the formulation of solving stationary axisymmetric spacetimes and focused on the 
example of doubly spinning BRs. The doubly spinning BRs are distinguished from the BRs with single 
spin by the new independent asymptotic conserved charge. 

This example shows that we can formulate the problem in a form suitable for numerical implemen- 
tation. As a first step, we have expanded the nonlinear equations for small second angular velocity 
and obtained the perturbative equations for the doubly spinning BRs. Separation of variables is not 
available for some set of equations, and we have integrated the perturbative equations numerically. 
The perturbations have been solved upto second order in small angular velocity. The constructed 
doubly spinning BRs are regular and free from the conical singularity. However, the absolute value of 
periodicity of axes, A(f> = Atp = 2ne, cannot be determined for each solution in the present approach. 
The ambiguity comes from the constant mode in e 2v ' . 

Although we have discussed the perturbations in the canonical coordinates, it is not necessary to 
work in these coordinates if we are interested only in the perturbative treatment of doubly spinning 
BRs. The corresponding analysis in the ring-coordinates will be more tractable than in the canonical 
coordinates. 
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To discuss the physical properties of doubly spinning BRs, we have evaluated the reduced entropy 
and reduced angular momentum for each perturbative solution. As peculiar to rotating black holes, the 
second spin associated to the S* 2 -sphere deceases the entropy for a fixed mass and angular momenta. On 
the other hand, the sum of reduced angular momenta decreases. Our construction of doubly spinning 
BRs is only for the branch of thin black rings, and it is unclear how the fat black rings behave with the 
second spin. An interesting question is whether the region where the Myers-Perry (MP) black holes 
and black rings overlap disappears or not. Presumably, the overlapped region will disappear for some 
range of second spin, because the branch of doubly spinning MP black holes with regular horizon 
leaves the overlapped region and its reduced entropy is bounded below (see Fig. [5|). If the doubly 
spinning BRs have smaller reduced entropy all the way down to its (expected) maximum spin, they 
can never reach the parameter region of MP black hole in some range of j 2 (or in three-dimensional 
space of {if, jf, ■§})■ 

The issue of dynamical stability of black rings and MP black holes are also important. The thin 
black rings will be unstable via Gregory-Laflamme instability along the S 1 direction [TJ [25], and the 
fat black rings have been though of more unstable \26\ 127] . Furthermore, the doubly spinning black 
rings presumably suffer from superradiant instabilities peculiar to a black hole with a rotating sphere 
[28]. 

We have also discussed interesting possible applications of numerics to constructing higher dimen- 
sional black rings and black rings in AdS. Although there are some possible technical subtleties involved 
in the approach, the direct approach to constructing new black rings is the great challenging. All the 
metric discussed above will be useful for our future study of black rings in higher dimensions and/or 
with the cosmological constant. We will be back these open problems in our future publications. 
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